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ABSTRACT 

We present an imaging study of the black hole binary candidate SDSS J1536+0441 (z « 0.3893), 
based on deep, high resolution VzK images collected at the ESO/VLT. The images clearly show an 
asymmetric elongation, indicating the presence of a companion source at ~ 1" (« 5 kpc projected 
distance) East from the quasar. The host galaxy of the quasar is marginally resolved. We find that 
the companion source is a luminous galaxy, the light profile of which suggests the presence of an 
unresolved, faint nucleus (either an obscured AGN or a compact stellar bulge). The study of the 
environment around the quasar indicates the occurrence of a significant over-density of galaxies with 
a redshift compatible with z w 0.4. This suggests that it resides in a moderately rich cluster of 
galaxies. 

Subject headings: quasars: individual (SDSS J153636. 22+044127.0) 



1. INTRODUCTION 

iBoroson fc Laueri (|2009f ) discovered peculiar features 
in the optical spectrum of SDSS J153636. 22+044127.0 
(hereon, SDSS J1536+0441), a z = 0.3893 radio quiet 
quasar. Three line systems are found at different red- 
shifts: a set of broad and narrow emission lines at z = 
0.3893, a set of only broad emission lines at z = 0.3727 
and a third set of very narrow abso rption lines at an 
interm ediate redshift (z = 0.3878). IBoroson fc Laueri 
(2009) interpreted the former two systems in terms of 
the broad line region emission around two massive black 
holes within the same galaxy, constituting a massive 
b lack hole binary ( BHB) w ith sub-parsec separation. 

IChornock et all (|2009allbh questioned the BHB sce- 
nario, ascribing the peculiar spectroscopic features of the 
quasar to a single active massive black hole with a partic- 
ular disk-like geometry of the broad line region, similar 
to Arp 102B. This scenario is also supported by the de- 
tection of a bump in the red wings of broad emission line 
profiles, at z w 0.4. 

Motivated by the c laim by IBoroson fc Laueri (2009), 
IWrobel fc Laon (120091) observed SDSS J1536+0441 us- 
ing the VLA, discovering two radio sources (VLA-A and 
VLA-B) separated by 0.97" within the quasar optical 
localization re gion. Independently, NIR observations led 
by our group (jDecarli et al.|[20"09bf ) revealed the occur- 
rence of a companion galaxy of the quasar, located in 
correspondence with VLA-B. This was subsequently con- 
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firmed in the optical wavelengths by ILauer fc Borosoiil 
(2009) using HST 

IWrobel fc Laon ((2009) and IDecarli et ai] l|2009b| ) sug- 
gested that each nucleus of the two galaxies is related 
to a single set of broad emission lines. The validity of 
this picture can be tested with high angular resolution 
spectroscopy, in order to prove if the two sets of broad 
emission lin es are spatially s e parated . ILauer fc Borosonl 
(2009) and IChornock etall (|2009bD presented optical 
spectra of SDSS J1536+0441. Both groups suggest that 
the regions responsible for the emission of the two sets of 
broad lines are spatially coincident. Nevertheless, both 
observations were carried out under non-optimal seeing 
conditions (larger than the angular separation of the two 
objects along the slit), hindering a clear comprehension 
of the contribution of the active nucleus in VLA-B to 
the spectrum of SDSS J1536+0441. 

The scenario of a sup erposition of two AG Ns has 
already been invoked by [Heckman et al] (|2009f ) to ex- 
plain the spectral properti es of another BHB candidate 
(SPSS J09 27+294 3: see IKomossa. Zhou fc Lul 120081: 
i Dotti et ail 120091 : iBogdanovic. Eracleous fc Sigurdssonl 
I2009D. but was sub s equent ly questioned by 
IDecarli. Reynolds fc DottH (|2009af) . IBoroson fc Laueri 
( 2009) noted that the expected number of chance quasar 
pairs with projected separations < 1" in the whole SDSS 
DR7 sample is only 0.3%. On the other hand, if the two 
objects are physically linked (e.g., if they lie in a rich 
galaxy cluster), such a number rises significantly (see 
the discussion in IWrobel fc Laorll2009f ). 

In order to study the properties of this peculiar object, 
and to probe the presence of a rich galaxy cluster in its 
environment, we performed deep, high-resolution imag- 
ing in the optical and near infrared wavelengths in an 
ESO Director Discretional Time programme (proposal 
ID: 282.B-5059). 

Throughout the paper, we adopt a concordance cos- 
mology with H a = 70 km/s/Mpc, Q m = 0.3, Q A = 0.7. 

2. OBSERVATIONS, DATA REDUCTION AND ANALYSIS 

K-band observations of the field around SDSS 
J1536+0441 were collected with the High Acuity Wide 
field K-band Imager (HAWK-I) at the ESO/ Very Large 
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Fig. 1.— K-band images of SDSS J1536+0441. Each panel is 
7" X 6" wide. North is up, East to the left, a) The quasar image, 
as observed in the K-band. The companion galaxy corresponding 
to VLA— B makes the quasar image clearly oblate towards the left . 
The contour plots of the radio detection by Wrobcl & Laor (2009) 
are also shown in the inset, plotted with the same size scale, b) 
Image of the companion source, after the subtraction of the quasar 
model, c) Image of the quasar, after the subtraction of the model 
of source B. d) Residuals after the subtraction of both the models. 
No obvious asymmetry in the residual map of the two sources is 
reported, suggesting that the two galaxies are fairly axisymmetric. 

Telescope (VLT) on April, 30, 2009. The total integra- 
tion time amounts to 48 minutes. HAWK-I is equipped 
with 4 CCDs of 2048x2048 pixels and a pixel-scale of 
0.1064"/pxl. Each CCD image 3.2 x 3.2 square arcmin 
(after trimming and the removal of partly observed re- 
gions). Usual jitter procedure was applied in order to col- 
lect broad band images in the NIR. Images were reduced 
with an IRAF 7 -based pipeline developed and tested in 
other programmes (e.g.. lKotilainen et al.f 2009). The see- 
ing was 0.75-0.80" (FWHM) throughout the whole inte- 
gration time. We estimate the Zero Point by comparing 
the instrumental magnitudes of field sources with the val- 
ues reported in the 2MASS database. The accuracy in 
the photometric calibration is 0.03 — 0.05 mag. Hereafter, 
we will focus on the Chip 2 image, the one overlapping 
onto the optical observations. 

V- and z-band images were collected with the Focal Re- 
ducer and low dispersion Spectrograph #2 (FORS2) at 
the ESO/VLT on May, 22 and 24, 2009. FORS2 mounts 
2 CCDs of 2048x4096 pixels. Our observations were car- 
ried out with the High-Resolution collimator with lxl 
binning, yielding a pixel-scale of 0.0632"/pxl. Offsets 
within a 10" wide box were applied in order to over- 
come the 3.5" gap between the two FORS2 CCDs. The 
final frames are 3.6 x 3.3 square arcmin wide. In the V- 
band, ten, 90 sec-long images were acquired, for a total 
integration time of 15 minutes. In the z-band, several 
120-sec long exposures were secured, yielding a total in- 
tegration time of 52 minutes. The seeing was excellent 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 



(FWHM=0.7" in z, 0.4" in V). The Zero Point of the z- 
band observation was estimated from the z magnitudes 
of field sources in the SDSS photometric catalogue. Sim- 
ilarly, for the V-band, we referred to the V magnitude 
estimates as derived from the SPSS q and r bands, fo l- 
lowing the conversion recipes by I Windhorst et all (|1991| ). 
The internal accuracy is around 0.04 mag. 

A composite image of our data in the 
three observed bands is available on-line at 
www . df m . uninsubria . it/astro/ qso_host/ql536/. 

The Ast ronomical Image Decomp osition and Analy- 
sis (AIDA, lUslenghi fc Falomoll2008f l software is used to 
model the Point Spread Functio n (PSF). The employed 
techn ique is widely discussed in iKotilainen et alj (|2007l 
2009), and it is briefly summarized here. The images of 
field stars are modeled with the superposition of 4 bidi- 
mensional gaussians mimicking the core plus an exponen- 
tial feature for the wings of the PSF. No significant PSF 
variation through the frame is reported in the K-band 
frame. On the other hand, FORS2 images are affected 
by severe distortions of the PSF shape throughout the 
image, making the PSF model less secure. 

3. RESULTS 
3.1. Decomposition of the quasar image 

Figure Q] shows that the image of SDSS J1536+0441 is 
clearly elongated East-wards, revealing the presence of a 
blended source (hereafter, source B) close to the quasar 
and coincident with VLA-B. In a comparison between 
the PSF model and the observed light profile of the two 
sources, it is apparent that they are both resolved (see 
Figure [2]). As quasar A is blended with its faint compan- 
ion, in order to study the presence of extended emission 
in both of them we adopt an iterative procedure to model 
the two sources: 1) First, we mask the Eastern side of the 
bright quasar (where the contribution from the compan- 
ion source may be relevant) and perform a 2D modelling 
with a nuclear+host galaxy profile convolved to our PSF 
model. 2) The model is subtracted from the observed 
frame, the residuals are masked out and we model the 
emission from the companion source. 3) The companion 
model is subtracted from the original observed frame, 
and we re- model the emission from the quasar. 4) Fi- 
nally, we re-perform the second step with the new model 
of the main source. 

In order to evaluate the uncertainties in the magni- 
tudes of nuclear and extended emissions, we scaled the 
Nuclear-to-Host luminosity ratio for a large set of models, 
and conservatively adopted as reference those bracketing 
the observed light profiles. Table [T] summarizes the re- 
sults of the fit procedure. The host galaxy of quasar A 
is significantly detected only in the K-band. On the con- 
trary, source B is well resolved. In all the bands, adding 
an unresolved nuclear component significantly improves 
the quality of the fit. As the unresolved-to-resolved lu- 
minosity ratio in the z band (~ rest-frame R band, see 
below) is around 0.3, we cannot assess whether it is a 
low-luminosity active nucleus (as suggested by the radio 
detection) or a stellar bulge. 

The observed V, z, K bands correspond to the rest- 
frame B, R and H. We compute the absolute magni- 
tudes in these bands as M x = m y — DM — A y + C x - y , 
where x refers to the rest- frame bands (B,R,H), y the 
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TABLE 1 

Best fit parameters in the modelling of the image of SDSS J1536+0441. (1) Band of the observations. (2) Total apparent magnitude of 
the model. (3) Apparent magnitude of the unresolved component. (4) Apparent magnitude of the host galaxy model. (5) Comparison 
between the x 2 obtained assuming a pure PSF and a galaxy+nucleus model. (6) The same as 5, comparing the cases without and with a 
nuclear component. (7) Rest frame band. (8) Absolute magnitude of the nucleus in the galaxy+nucleus decomposition. (9) The same as 

8, for the galaxy. 
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Fig. 2.— The K-band light profiles of SDSS J1536±0441 (A) and 
its companion source (B) as compared to the best fit model (solid 
line), the PSF model (dotted line) and the galaxy model (dashed 
line). The light profile of quasar A shows deviations from the model 
at large r adii, suggesting that the host galaxy may present stellar 
shells (cfr Canalizo et ah 2 0071) . 



observed ones (V,z,K), DM is the distance modulus, 
A y is the Galactic extinction and C x —y is a term ac- 
counting for filter and fc-correction. In our cosmologi- 
cal framework, z = 0.3893 yields DM = Aim. The 
Galactic extinction is derived from the HI maps by 
ISchlegel. Finkbeiner fe Davisl (|1998l) . T he fc-correction 
comes from the galaxy templates by iMannucci et al 



2001) and the quasar template by iVanden Berk et al 



2001). The resulting magnitudes are reported in table[TJ 



Source A has host galaxy absolute magnitudes consistent 
with the typical quas ar host galaxy luminosities (e.g., 
iKotilainen et all 120091 ) and it is ~ 2 mag brighter than 
the char acteristic luminosity of quiescent galaxies at this 
redshift (jCirasuolo et al.ll200l . The companion source 
is a galaxy ~ 1 mag fainter (cfr also the M g = —21.4 es- 



timate reported bv lLauer fc Borosorj|l2009T ) with a much 
fainter, red nucleus. 

We note that, if the B-band flux of the unresolved 
source is interpreted in terms of an active galactic nu- 
cleus, source B should be classified as r adio loud, accord- 
ing t o the standard definition (e.g., iKellermann et ail 
1989) based on the radio-to-optical flux ratio (here, R ~ 
150, assuming the radio flux found by iWrobel fe Laorl 
[20M ). 

3.2. The properties of the environment 

In order to investigate the clustering properties of the 
environment of this quasar, we search for galaxies in its 
field. We consider the 3.2 x 3.2 square arcmin (« 1 x 1 
Mpc at z = 0.4) region that has been gathered in all 
the bands. We produce an inventory of sources ex- 
ceeding a 5-er detection threshold with respect to the 
background count rms t hrough the software SExtractor 
(|Bertin fc Arnoutill996| ). The corresponding limit mag- 
nitudes are rriy K = 25.7, 23.1 and 19.8. We then match 
the compilations of sources detected in the various bands. 
The complete list consists of 200 sources detected in the 
V-band, 165 in z (153 of which also appearing in the V- 
band catalogue) and 87 in K (80 of which also detected 
in V). A total of 69 sources are detected in all the bands. 

Among them, we select galax y candidates using the 
color-color diagram proposed bv lBlanc et alj (|2008l ) and 
adapted to our photometric system (see Figure [3]). We 
assume that a galaxy has: 



(z-K) > 0.3(V- z) + 1.85 



(1) 



Out of 69 sources, 34 match this color criterion. We note 
that, given our sensitivity limits, we expect to detect 
~ 30 stars in our frame in all the bands at this Galactic 
latitude 8 . Similarly, if we conside r the number density 
of general field ga laxies (e.g., from lCirasuolo et"al~ll2008l : 
iBlanc et al.l 120081) . 14 candidates brighter than our sen- 
sitivity limits are expected. This indicates that the field 
around t he quasar is over-dense . From the galaxy tem- 
plates by IMannucci et all (|2001[ ). we infer that a z « 0.4 

8 Based on the Galaxy simulations per- 
formed with the software TRILEGAL avail able at 
http : //stev . oapd . inaf . it/cgi-bin/trilegal HGirardi et all 
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Fig. 3. — The VzK color diagram for the objects in our anal- 
ysis. Symbol size maps the object brightness. 'Ell', 'Sc' and 
'QSO' la bel the expectation from the elliptical and Sc galaxy tem- 
plates by Marmucci ct al. (2001) and from the quasar template by 
IVanden Berk et alj 11200 ll) . redshifted to z = 0.4. The solid lines 
show how galaxy colors change from z = to z = 0.4. The typical 
error box is shown as a cross in the corner of the frame. The color 
cut in equation[T]is shown as a dashed line. A number of sources lie 
within the expectation for galaxies at z ~ 0.4 (the vertical, dotted 
lines), suggesting the presence of a galaxy cluster. 

galaxy has V-z ra nging between 1.5 and 2.5 (see also 
ICsabai eT al. 2003); thus, 18 out of 34 galaxy candidates 
are consistent with being at z ~ 0.4. The correspond- 
ing absolute magnitudes of these objects, assuming that 
they are galaxies at the same redshift as the quasar, are 
(M B > = -22.8, (M R ) = -20.5 and (M H ) = -19.2 (with 
rms ~ 1 ma g), consistent with ty pical fairly bright galax- 
ies (cfr, e.g. IGavazzi et al. 2003). It is interesting to note 
that, out of 11 well-resolved sources in the K-band, 10 
are classified as galaxies in our color selection, and 9 of 
them are consistent with having z m 0.4. 

We compute the number density of these galaxies as a 
function of the angular distance from the quasar in order 
to evaluate the occurrence of any clustering (see Figure 
13]) ■ We find a significant excess of the galaxy density 
within < 200 kpc from SDSS J1536+0441 with respect 
to average field. Thus the quasar is found in a moderately 
rich cluster of galaxies. 

4. DISCUSSION 

In this Letter, we presented a high-resolution, deep 
multi-color imaging study of the alleged binary quasar 
SDSS J1536+0441. In addition to the extended emission 
of the host galaxy, we detect a bright companion galaxy 
at ~ 1" from the quasar. The quasar and the compan- 
ion are located in c orrespondence w i th two distinct radio 
sources observed bv lWrobel fc Laorl (|2009| ). showing that 
both galaxies host an AGN. 

Our deep images unveil the presence of a rich galaxy 
cluster. The existence of such a massive structure solves 
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Fig. 4. — The number density of the galaxy candidates consistent 
with z = 0.4 as a function of the projected distance from SDSS 
J1536+0441. For a comparison, the expected number density of 
background galaxies is plotted as a dashed line. A clear excess of 
galaxy candidates within ~ 200 kpc from the quasar is observed, 
supporting the cluster scenario. 

the apparent inconsistency between the observation of 
two radio sources within < 1" and the expected num- 
be r of such pairs in the S PSS catalogue, as di s cusse d 
in iBoroson fe Lauerl ()2009l ) and IWrobel fc Laorl (2009). 
Furthermore, the potential well of such a rich galaxy clus- 
ter could in principle explain the shift in redshift between 
the two sets of emission lines. This consideration favours 
a superposition picture to explain the spectrum of SDSS 
J1536+0441. 

However, the red color and the faintness of the 
nuclear optical continuum observed in the compan- 
ion galaxy suggest that its AGN emission is obscured. 
As a consequence, the AGN of the companion galaxy 
should not contribute much to the broad emission lines 
present in the SDSS spe ctrum, in agreement wit h the 
indications discu s sed in iChornock et alJ (|2009bl ) and 
lLauer fc Borosonl (|2009l ). In this case, the two sets of 
broad emission lines should be emitted by the quasar, 
consistently with the BHB or the double peaked emitter 
scenarios. The external regions of the companion galaxy 
could also naturally explain the narrow absorption lines 
observed in the spe ctrum of SDSS J 1 536+0 441, as al- 
ready suggested by lLauer fc Boroson (2009). The ori- 
gin of_Jhe_j^d_j2ump_Jn_the broad line profiles reported 
bv lChornockeTaTI (|2009bh remains unclear both in the 
BHB and the superposition scenarios. 

The effective role of the companion in the spectrum of 
SDSS J1536+0441 and of the surrounding galaxy cluster 
is yet to be clarified. The superposition model, where the 
companion emits one of the two sets of emission lines is 
in possible conflict with the presently available spectro- 
scopic data. However, the discovery of a possible nu- 
cleus in source B, ma tched with its radio detection by 
IWrobel fc Laorl (|2009f ). may support this picture. Cru- 



Probing the nature of SDSS J1536+0441 



5 



cial tests for this interpretation can be derived: from 
UV observations of hig h-ionization lines (as discussed in 
iChornock et al.ll2009bf) . high angular resolution optical 
spectroscopy and high-resolution X-ray images. 
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